: 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

COMPUTER  AIDED  DATA  ACQUISITION  AND  CONTROL 

OF  AN   INTERNAL   COMBUSTION  ENGINE 

by 

Bryan  R.   Oakes 

March  1984 

Thesis 

;  Advisor:                                                    W.G. 

Culbreth 

Approved  for  public  release;  distribution  unlimited 


T215658 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whit  Diti  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.    REPORT  NUMBER 


2.  GOVT  ACCESSION  NO 


3.     RECIPIENT'S  CATALOG  NUMBER 


4.    TITLE  (and  Subtitle) 


5.     TYPE  OF   REPORT  A   PERIOD  COVERED 


Computer  Aided  Data  Acquisition  and  Control    of 
an   Internal    Combustion  Engine 


Master's  Thesis; 
March   1984 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.     AUTHORS 

Bryan  R.  Oakes 


8.  CONTRACT  OR  GRANT  NUMBERf*,) 


••  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Naval  Postgraduate  School 
Monterey,  California  93943 


10.     PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  4   WORK  UNIT  NUMBERS 


II.    CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Naval  Postgraduate  School 
Monterey,  California  93943 


12.     REPORT   DATE 

March  1984 


>3.     NUMBER  OF  PAGES 


14.    MONITORING  AGENCY  NAME  S   AOORESSfM  different  from  Controlling  Office) 


106    pages 


15.     SECURITY  CLASS,   (of  thla  report) 

Unclassified 


15*.     DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


l«.    DISTRIBUTION  STATEMENT  (of  thle  Report) 

Approved  for  public  release;  distribution  unlimited 


17.    DISTRIBUTION  STATEMENT  (of  the  ebetrect  entered  In  Block  20,  It  dltlerent  from  Report) 


It.    SUPPLEMENTARY  NOTES 


It.    KEY  WORDS  (Continue  on  reeeree  aid*  If  nocooemry  and  Identity  by  block  number) 

Data  Acquisition,  Digital  Control,  Computer  Control,  Microprocessor, 
Engine  Control 


l$-     ABSTRACT  (Continue  on  reveree  elda  (I  neeaeeary  and  Identity  by  block  number)         .,_.-—«  ,  _ 

Two  microcomputers,   a  Digital    Equipment  Corporation  VT-103  and  an  Octagon 
Systems  Corporation  SYS-1,  were  interfaced  with  a  General    Motors  Model    1-53X3 
single  cylinder  diesel   engine  to  provide  automated  data  acquisition  and  data 
reduction  and  engine  control   while  running  classroom  experiments.     Engine  inlet 
and  exhaust  temperatures  and  pressures,   RPM,  torque  and  fuel    flow  along  with 
additional    engine  parameters  can  be  measured  with  the  aid  of  a  computer^ 
Values  for  parameters  such  as  specific  fuel    consumption,   air-to-fuel    ratio, 
and  thermal   efficiency  can  then  be  computed  and  both  measured  and  computed 


do,; 


'ST*  1473 


EDITION  OF  I  NOV  «S  IS  OBSOLETE 
S/N  0102- LF- 014-6601 


I     SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Datm  Entered' 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Dalm  Enrarad) 


20.  (Continued) 

values  can  be  displayed  and  recorded.  The  Digital  Equipment  Corporation 
computer  is  the  parent  computer  and  is  used  interactively  for  data 
acquisition  and  data  reduction  and  for  feedback  control  through  the 
Octagon  computer  to  which  it  is  linked.  The  Octagon  computer  is  used 
exclusively  for  control  of  the  engine  speed  and  load. 


$<N  0)02-  LF-  0)4-6601 


2      SECURITY  CLASSIFICATION  OF  THIS  P AG€.(Wh»n  Dmtm  Knfrmd) 


Approved  for  public  release;  distribution  unlimited 


Computer  Aided  Data  Acquisition  and  Control  of  an 
Internal  Combustion  Engine 


by 


Bryan  R.  Oakes 

Lieutenant,  United  States  Navy 

B.A.,  Linfield  College,  Oregon,  1975 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
March  1984 


DUD 
ABSTRACT 


• 


Two  microcomputers,  a  Digital  Equipment  Corporation 
VT-103  and  an  Octagon  Systems  Corporation  SYS-1 ,  were 
interfaced  with  a  General  Motors  Model  1-53X3  single  cylin- 
der diesel  engine  to  provide  automated  data  acquisition  and 
data  reduction  and  engine  control  while  running  classroom 
experiments.   Engine  inlet  and  exhaust  temperatures  and 
pressures,  RPM,  torque  and  fuel  flow  along  with  additional 
engine  parameters  can  be  measured  with  the  aid  of  a 
computer.   Values  for  parameters  such  as  specific  fuel 
consumption,  air-to-fuel  ratio,  and  thermal  efficiency  can 
then  be  computed  and  both  measured  and  computed  values  can 
be  displayed  and  recorded.   The  Digital  Equipment 
Corporation  computer  is  the  parent  computer  and  is  used 
interactively  for  data  acquisition  and  data  reduction  and 
for  feedback  control  through  the  Octagon  computer  to  which 
it  is  linked.   The  Octagon  computer  is  used  exclusively  for 
control  of  the  engine  speed  and  load. 
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I.   INTRODUCTION 

A  computer  aided  data  acquisition,  data  reduction  and 
engine  control  system  was  constructed  and  tested  on  a  diesel 
engine  test  stand  used  for  laboratory  experimentation.   The 
engine  used  was  a  General  Motors  model  1-54X3  diesel  engine. 
The  engine  was  coupled  to  a  DC  motor  generator  which  was 
used  as  a  dynamometer  to  place  a  variable  load  on  the 
engine.   Power  generated  was  dissipated  in  a  resistor  load 
bank . 

In  the  past,  whenever  the  unit  was  used  for  classroom 
experimentation,  the  recording  of  engine  parameters  and  the 
reduction  of  the  data  was  done  by  hand.   A  technician  was 
also  required  to  run  the  engine  during  testing.   In  order  to 
more  efficiently  and  accurately  record  data  and  to  allow  a 
means  of  immediately  reducing  and  displaying  the  engine's 
performance  parameters,  a  microprocessor-based  data 
acquisition  and  data  reduction  system  was  designed.   In 
addition,  in  order  to  reduce  the  need  for  a  technician  and 
to  allow  the  engine  to  be  brought  to  a  desired  operating 
point  for  study  more  quickly  and  easily,  a  microprocessor- 
based  control  system  was  designed  for  both  the  engine  and 
dynamometer . 

An  Octagon  SYS-1  8-bit  microcomputer  was  chosen  as  the 
control  computer.   Digital-to-analog  and  analog-to-digital 
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converters  were  built  to  allow  the  computer  to  interface 
with  transducers  attached  to  the  engine/dynamometer  unit. 
The  transfer  functions  for  both  the  diesel  engine  and  the 
dynamometer  were  experimentally  determined  and  used  to 
design  an  acceptable  control  algorithm  for  use  by  the 
control  computer. 

The  use  of  the  microcomputer  in  the  feedback  control 
loops  allowed  the  use  of  almost  any  control  action.   A 
proportional-plus-integral-plus-differential  (PID)  control 
action  was  chosen  [Ref.  11  [Ref.  2].   The  use  of  a  PID 
control  action  allows  the  designer  of  a  control  system  a 
measure  of  control  over  the  speed  of  response  of  the  system 
to  a  step  input  and  gives  acceptable  steady  state  accuracy 
when  used  with  the  hardware  available.   The  values  for  the 
three  constants  associated  with  each  PID  controller;  engine 
and  dynamometer;  were  initially  determined  using  Bodie 
diagrams  of  the  open  loop  frequency  response  curves.   A 
short  control  program  was  written  for  the  control  computer 
using  the  analytically  determined  values  for  the  constants 
and  the  program  was  tested  after  being  installed  in  the 
system.   The  final  values  for  the  constants  used  in  the 
control  program  were  experimentally  determined. 

For  data  acquisition  and  data  reduction,  a  Digital 
Equipment  Corporation  VT-103  desktop  computer  was  used.   The 
VT-103  was  interfaced  with  the  SYS-1  computer  to  allow  user 
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input  of  desired  engine  RPM  and  load.   The  Digital  computer 
included  a  16  channel  analog-to-digital  converter  that 
provided  a  means  of  linking  the  computer  to  the  sensing 
devices  installed  on  the  diesel  engine  test  stand.   An 
interactive  FORTRAN  program  was  written  that  allowed  the 
user  to  run  a  variety  of  experiments  on  the  system  while 
controlling  the  engine  RPM  and  load. 
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II.   HISTORY 

The  use  of  microprocessors  in  control  and  data 
acquisition  applications  has  grown  dramatically  since  the 
early  1970's.   Computers  were  first  used  for  control  in  the 
preceeding  decade.   Early  in  the  60's,  after  computers  began 
to  be  used  in  business  applications,  the  advantages  of 
computers  as  a  process  control  device  first  began  to  become 
apparent . 

Computers  had  several  advantages  over  conventional 
control  devices.   Until  the  advent  of  computers,  process 
control  was  generally  broken  down  into  as  many  simple 
control  systems  as  possible.   The  only  device  capable  of 
reading  a  large  number  of  parameters  and  making  control 
decisions  based  on  all  of  them  was  a  human.   Computers,  both 
analog  and  digital,  because  they  could  perform  complex 
computations  with  constant  accuracy  and  at  high  speed, 
became  desirable  as  integral  parts  of  complex  manufacturing 
and  industrial  process  control  systems.   In  most  cases,  even 
though  they  were  expensive,  they  were  more  cost  efficient 
and  accurate  than  the  human  operators  they  replaced. 

During  this  same  time  period  the  advantages  of  digital 
over  analog  computers  became  apparent.  Almost  any  desired 
increase  of  accuracy  could  be  achieved  using  digital 
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computers  with  relatively  little  increase  in  cost  due  to  the 
increased  accuracy.   It  was  almost  impossible  to  increase 
the  accuracy  of  analog  computers  without  greatly  increasing 
the  cost  of  the  computer. 

Early  in  the  1970's  the  first  microcomputer  chips  became 
available  on  the  market.   INTEL  marketed  the  first  8-bit 
microprocessor  in  1972.   The  first  applications  of  the  new 
microprocessors  were  in  hand  held  calculators,  but 
inexpensive  microcomputers  (8-bit)  and  minicomputers 
(16-bit)  based  on  the  ne.w  chips  soon  followed. 

Since  the  early  70's,  the  applications  of  micro- 
processors has  grown  rapidly.   The  cost  of  the  computers 
dropped  exponentially  during  the  decade  and  at  the  same  time 
reliability  improved.   Microcomputers  are  now  used  in  such 
diverse  control  areas  as  chemical  process  control,  aircraft 
control  and  machining.   In  all  of  these  applications  the 
microprocessor  samples  parameters  and  makes  control 
decisions  based  on  the  samples. 

The  first  applications  of  computer  process  control 
required  new  methods  of  analysis  of  the  control  problem. 
Digital  computers  were  unable  to  sample  control  parameters 
continuously,  instead,  the  computer  would  sample  once  during 
a  discrete  time  period.   This  led  to  the  development  of  new 
methods  for  analyzing  the  control  problem.   Continuous 
methods  of  analysis  can  still  be  used  in  some  cases  but  they 
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may  not  be  as  accurate  as  some  of  the  discrete  time  period 
methods  [Ref .  3  ]  • 

One  of  the  first  applications  of  computers  was  in  data 
reduction.   When  first  employed  in  this  capacity,  data  that 
had  been  recorded  by  hand  was  entered  into  a  computer  which 
would  then  perform  the  necessary  calculations.   The  tech- 
nology that  made  computer  process  control  possible  also  made 
computer-aided  data  acquisition  possible.   The  interfaces 
between  the  analog  real  world  and  the  digital  computers 
necessary  for  the  measurement  of  control  parameters,  also 
allowed  computers  to  record  data  in  addition  to  reducing 
data  . 
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III.   DESCRIPTION  OF  SYSTEM  HARDWARE 

The  design  of  both  the  data  acquisition/data  reduction 
subsystem  and  the  control  subsystem  evolved  around  existing 
hardware.   An  effort  was  also  made  to  use  existing  hardware 
when  available  and  to  limit  the  cost  of  any  item  that  needed 
to  be  purchased  or  manufactured.   The  diesel  engine  and 
dynamometer  have  been  in  use  since  1960.   The  system  was 
upgraded  in  1979  to  include  more  up-to-date  transducers  and 
measurement  devices  and  digital  readouts  of  many  of  the 
student  recorded  parameters  and  a  new  control  console  for 
use  by  the  technician  running  the  engine.   At  the  time  of 
the  update,  provisions  were  made  to  install  a  control  system 
and  a  data  acquisition  system  when  it  became  available.   The 
layout  of  the  engine  and  dynamometer  and  associated 
equipment  is  shown  in  Figure  1. 

A.   ENGINE  AND  ASSOCIATED  HARDWARE 

The  engine  used  in  the  experiment  was  a  General  Motors 
Model  1-53X3  two  stroke-cycle  diesel  engine.   The  engine  was 
coupled  to  a  DC  motor  generator  that  could  be  used  as  a 
dynamometer  by  feeding  the  generated  current  to  a  load  bank 
to  be  dissipated  as  heat  or  that  could  be  used  as  an 
electric  motor  to  start  the  diesel. 

A  total  of  eleven  parameters  associated  with  the  diesel 
engine  needed  to  be  measurable  by  the  data  acquisition 
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computer.   This  required  that  the  output  from  all  trans- 
ducers be  converted  to  a  voltage  range  that  was  compatible 
with  transistor-transistor  logic  (TTL)  devices.   Many  of  the 
existing  transducers  had  signals  already  conditioned  to  be 
used  by  the  computer . 

Six  Newport  pyrometers  were  installed  to  measure  intake 
and  exhaust  gas  temperatures,  coolant  temperatures,  fuel 
temperature  and  ambient  air  temperature.   The  pyrometers  had 
a  signal  output  that  varied  between  0  and  +5  volts  DC  and 
required  no  additional  electronics  to  be  used  with  the  data 
acquisition  computer. 

Pressures  had  in  the  past  been  measured  using  mercury 
and  water  manometers.   Three  strain  gage  differential 
pressure  transducers  were  purchased  which  included  the 
electronic  circuitry  necessary  to  amplify  the  strain  gage 
output  to  the  range  of  0  to  +5  volts. 

Another  differential  pressure  transducer  was  purchased  to 
measure  torque.   The  existing  torque  measuring  device  con- 
sisted of  a  bellows  attached  to  a  lever  arm  linked  to  the 
dynamometer.   A  pressure  sensing  Bourdon  tube  analog  meter 
converted  the  pressure  from  the  bellows  to  a  torque  reading. 
The  differential  pressure  transducer  was  installed  in  series 
with  the  existing  gage  and  calibrated  to  measure  lbf-ft. 
Engine  RPM  was  measured  with  a  magnetic  pickup  on  the  dynamo- 
meter shaft.   Fuel  flow  was  measured  with  a  turbine  flow 
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meter.   Both  of  these  transducers  produce  variable  frequency 
signals.   Instrumentation  amplifiers  were  used  to  convert  the 
frequency  output  to  a  variable  voltage  for  use  by  the  computer 
Provisions  were  also  made  to  measure  coolant  mass  flow  but  a 
transducer  was  yet  to  be  installed.   A  list  of  all  of  the 
sensors  and  their  characteristics  is  found  in  Appendix  A. 

B.   DATA  ACQUISITION  AND  DATA  REDUCTION  SYSTEM  HARDWARE 

A  Digital  Equipment  Corporation  VT-103  microcomputer  was 
available  to  be  used  for  data  acquisition/reduction,  or 
engine  control,  or  both.   The  VT-103  was  a  16-bit  micro- 
computer based  on  DEC'S  LSI -11  microprocessor.   The  computer 
came  with  an  8  slot  bus  that  could  be  configured  with  a 
number  of  hardware  options  dependent  on  the  needs  of  the 
user.   It  could  be  programmed  in  FORTRAN  and  had  an 
operating  system  with  a  number  of  subroutines  callable  from 
FORTRAN  programs  that  made  it  very  easy  to  access  various 
hardware  options  without  the  need  of  software  drivers; 
analog-to-digital  converters,  digital-to-analog  converters, 
serial  interfaces  and  parallel  interfaces. 

The  decision  was  made  to  dedicate  the  VT-103  to  data 
acquisition  and  data  reduction.   An  interface  was  eventually 
provided  to  connect  the  control  computer  with  the  VT-103. 
To  properly  run  a  variety  of  experiments  and  to  allow 
the  user  a  certain  amount  of  flexibility,  it  was  anticipated 
that  an  interactive  data  acquisition/reduction  program  would 
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be  needed.   The  flexibility  and  ease  of  programming  made  the 
VT-103  ideal  for  this  application.   It  was  originally  felt 
that  the  control  computer  would  also  need  an  interactive 
program  to  allow  the  user  a  means  of  changing  the  operating 
parameters  of  the  engine.   Instead,  the  two  computers  were 
interfaced  through  an  asynchronous  parallel  line.   This 
greatly  simplified  the  programming  required  by  the  control 
computer  and  allowed  all  of  the  control  functions  in 
addition  to  all  of  the  data  acquisition  functions  to  be 
performed  by  the  operator  on  one  terminal. 

In  the  final  configuration  the  VT-103  computer  required 
a  means  of  storing  the  final  program,  a  terminal  for 
interactive  display,  a  printer,  a  multi-channel  analog-to- 
digital  converter,  and  a  parallel  port  for  interfacing  with 
the  control  computer.   Figure  2  show  the  configuration  of 
the  VT-103  backplane  with  all  of  the  required  options. 

The  LSI-11  microprocessor  with  4K  by  16-bit  random 
access  memory  (RAM)  resided  in  slots  1  and  2  of  the 
backplane.   Slot  3  held  a  multifunction  module.   This  module 
contained  an  additional  4K  by  16-bit  read/write  memory  and 
two  asynchronous  line  interfaces.   The  interfaces  were  used 
to  connect  the  microprocessor  to  the  terminal  and  to  the 
internal  tape  drives.   Slot  4  contained  an  additional  32K  by 
16-bit  RAM.   The  serial  port  for  the  printer  resided  in  slot 
5;  and  the  parallel  port  for  the  control  computer  resided  in 
slot  6.   Slots  7  and  8  were  for  the  16  channel  by  12-bit 
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analog-to-digital  (A/D)  converter.   Complete  descriptions  of 
the  boards  and  their  applications  can  be  found  in  Digital's 
microcomputer  handbooks  [Ref.  4]  and  [Ref.  53. 

C.   CONTROL  SYSTEM  HARDWARE 

The  hardware  required  for  the  control  system  did  not 
need  to  be  nearly  as  complex  as  the  hardware  required  for 
data  acquisition/reduction.   For  the  engine  and  dynamometer 
control,  only  two  parameters  needed  to  be  measured;  engine 
RPM  and  torque;  and  two  control  voltages  needed  to  be 
outputted;  one  to  the  diesel  engine  throttle  servo  and  one 
to  the  dynamometer  field  current  controller.   Both  the  servo 
and  field  current  controller  were  already  configured  to  be 
actuated  by  voltages  in  the  range  of  TTL  devices.   Although 
a  terminal  was  not  required  in  the  final  configuration,  one 
was  installed  while  programming  the  control  computer  and 
while  making  dynamic  adjustments  to  the  program. 

An  Octagon  SYS-1  microcomputer  was  available  for  use 
as  the  control  computer.   The  SYS-1  was  based  on  the  INS8073 
microcomputer  chip.   It  was  an  8-bit  computer  that  came  with 
an  RS232  input/output  port  and  a  three  channel  peripheral 
interface  chip  (PIO).   In  addition,  all  of  the  INS8073  pins 
were  available  at  a  44  pin  edge  connector  that  functioned  as 
a  bus.   All  of  the  circuitry  for  the  SYS-1  was  located  on 
one  5x7  circuit  board.   A  schematic  of  the  computer  can  be 
found  in  the  SYS-1  User's  Guide  [Ref.  6]  or  in  an  article  by 
McKown  and  Sarns  [Ref.  73 •   In  addition  to  the  microprocessor 
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the  board  contained  4K  of  RAM,  a  utility  program  held  in  a 
1K  byte  eraseable  programmable  memory  (EPROM)  chip  and 
another  EPROM  that  could  be  programmed  by  the  user. 

The  SYS-1  did  not  come  with  an  A/D  converter  or  a 
digital-to-analog  (D/A)  converter.   Both  were  needed  in  the 
final  system.   The  A/D  converter  was  necessary  to  read  the 
output  from  the  RPM  and  torque  measuring  transducers  and  the 
D/A  converters  to  convert  the  digital  control  signals  to 
control  voltages.   Rather  than  purchase  the  necessary 
hardware  it  was  felt  that  they  could  as  easily  be  built. 
Figure  3  is  the  schematic  for  the  D/A  converter  that  was 
designed  to  be  used  with  the  SYS-1 .   Figure  4  is  the  sche- 
matic of  the  A/D  converter.   Each  was  built  on  a  separate 
circuit  board  and,  except  for  the  voltage  inputs  to  the  A/D 
board  and  the  voltage  outputs  from  the  D/A  board,  used  the 
same  layout  for  the  44-pin  edge  connectors  as  the  SYS-1 . 

The  A/D  converter  was  an  8  channel  by  8-bit  converter. 
Only  two  channels  are  presently  used.   The  D/A  converter  had 
three  channels  and  was  also  an  8-bit  converter.   Only  two  of 
its  channels  were  used. 

The  three  channel  PID  was  used  to  interface  the  SYS-1 
with  the  data  acquisition  computer.   Two  of  the  8-bit 
channels  were  connected  to  the  16-bit  parallel  port  of  the 
VT-103.   The  8  high  bits  from  the  VT-103  parallel  port  were 
connected  to  channel  B  of  the  SYS-1  and  the  8  low  bits  were 
connected  to  channel  A. 
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IV.   DATA  ACQUISITION  AND  DATA  REDUCTION  THEORY 

The  data  acquisition  system  was  based  on  the  Digital 
VT-103  minicomputer  which  used  a  12-bit,  16-channel  analog- 
to-digital  (A/D)  converter  to  sample  various  parameters 
associated  with  engine  performance.   Prior  to  designing  the 
software  associated  with  the  system,  decisions  were  required 
on  sampling  rate  and  the  number  of  samples  and  an  analysis 
needed  to  be  made  to  determine  the  accuracy  of  the  data 
acquisition  system. 

A.   ANALOG-TO-DIGITAL  CONVERSION 

A  computer's  central  processing  unit  performs  all  of  its 
functions  using  binary  numbers.   In  general,  a  computer 
cannot  sense  or  use  any  analog  signal  directly.   The  analog- 
to-digital  converter  acts  as  the  interface  between  an  analog 
device  such  as  a  transducer  and  the  central  processing  unit 
of  a  computer.   The  analog-to-digital  converter  does  this  by 
converting  an  analog  signal  to  a  binary  number  that 
corresponds  to  the  magnitude  or  strength  of  the  signal.   An 
8-bit  analog-to-digital  converter  will  convert  an  analog 
voltage  to  a  binary  number  between  00000000  and  11111111 
(0  and  255  decimal).   A  12-bit  analog  converter  can  convert 
an  analog  voltage  into  an  integer  number  between  0  and  4095 
decimal . 
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The  VT-103  used  a  16  channel  12-bit  analog-to-digital 
converter.   Only  10  of  the  channels  were  required  to  measure 
the  diesel  engine  parameters.   All  of  the  parameters  were 
converted  to  voltages  between  0  and  +  5  volts  DC  by  the 
transducers  and  their  associated  signal  conditioning 
devices.   The  A/D  converter  used  in  this  system  could  accept 
signals  in  the  range  of  -5  to  +5  volts.   A  signal  of  -5 
volts  would  be  converted  to  decimal  0  while  a  signal  of  +  5 
volts  will  be  converted  to  decimal  4095.   It  should  be  noted 
that  if  the  signal  had  a  potential  of  0  volts  the  computer 
will  read  2048,  not  0. 

The  A/D  converter  used  with  the  VT-103  was  a  successive 
approximation  analog-to-digital  converter.   It  converted  an 
analog  signal  by  successively  comparing  the  signal  to 
various  reference  voltages  through  a  series  of  resistors. 
The  12-bit  converter  used  with  the  VT-103  required  twelve 
successive  comparisons  to  convert  one  signal  to  a  binary 
number . 

The  resolution  of  the  12-bit  converter  was  equal  to 
1/4095  times  the  acceptable  input  range  of  the  analog 
signal.   Since  the  VT-103  A/D  converter  could  accept  signals 
between  -5  and  +5  volts,  the  resolution  of  the  converter 
used  in  this  system  is  equal  to, 


(1/4095)  x  10  volts  =  +0.00244  volts 


(eqn.  4.1) 
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To  calculate  the  resolution  of  the  converter  in 
engineering  units,  the  resolution  of  the  converter  is 
divided  by  the  resolution  of  the  transducer  given  in 
volts/unit.   As  an  example,  the  resolution  of  the  Newport 
pyrometers  was  .001  volts/degree  F.   This  gives  a 
temperature  resolution  of  the  A/D  converter  of, 

.00244/. 001  =  +2.44  degrees  F         (eqn.  4.2) 

A  more  practical  method  of  computing  the  resolution  of  the 
A/D  converter  is  through  the  use  of  the  data  that  was 
collected  for  the  converter  calibration.   The  range  measured 
in  engineering  units  can  be  divided  by  the  range  of  A/D 
converter  values  to  give  the  same  result.   Using  the  same 
temperature  example,  temperatures  were  measured  between  60 
and  638  degrees  F,  a  range  of  578  degrees  F.   The  A/D 
converter  values  ranged  from  2068  to  2297,  a  range  of  229. 
Dividing  the  temperature  range  by  the  A/D  range  gives  a 
resolution  of  +2.52  degrees  F.   The  resolution  of  every 
channel  used  in  this  system  is  given  in  Appendix  B  along 
with  the  data  used  to  calibrate  the  analog-to-digital 
converter . 

B.   NYQUIST  CRITERION  AND  NUMBER  OF  SAMPLES 

Two  additional  factors  are  normally  considered  when 
designing  a  computer-aided  data  acquisition  system;  the 
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frequency  that  data  is  sampled  and  the  number  of  represen- 
tative samples  to  be  taken  and  averaged  for  any  set  of 
calculations . 

The  Nyquist  criterion  gives  a  good  indication  of  the 
frequency  that  data  should  be  sampled.   According  to  the 
criterion,  samples  should  be  taken  at  a  frequency  at  least 
twice  that  of  the  highest  natural  frequency  of  interest.   A 
good  explanation  of  the  need  for  Nyquist  criterion  can  be 
found  in  Pearson's  article  on  processing  analog  signals 
which  appeared  in  Popular  Electronics  [Ref.  8].   Sampling  at 
twice  the  highest  frequency  prevents  signal  aliasing; 
interpreting  the  signal  as  a  lower  frequency  signal  than  it 
really  is.   Since  all  data  for  the  diesel  is  taken  while  the 
diesel  is  running  at  a  constant  RPM  and  load,  i.e.,  steady 
state,  there  is  no  danger  of  sampling  the  data  too  slowly. 
In  actuality,  the  computer  samples  data  so  quickly  that 
problems  associated  with  low  frequency  sampling  do  not 
occur  . 

The  second  factor  considered  is  the  number  of  samples 
needed  to  be  taken  in  order  to  calculate  the  value  of  any 
parameter  with  reasonable  accuracy.   No  criterion  was  used 
to  determine  the  number  of  samples.   Instead,  the  program 
was  written  to  take  an  arbitrary  number  of  samples  (100)  and 
compute  the  average  value  and  the  standard  deviation  of  the 
100  samples  for  each  channel.   For  most  channels,  the 
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standard  deviation  was  consistently  less  than  1/4095  or 
0.02%   For  engine  RPM  the  standard  deviation  varied  slightly 
around  13  out  of  4096  or  0.24%.   This  accuracy  corresponded 
to  +12  RPM. 

C.   CALIBRATION  OF  THE  DIGITAL-TO-ANALOG  CONVERTER 

For  each  data  run,  after  the  computer  sampled  and 
averaged  the  data,  the  decimal  numbers  representing  the 
amplitude  of  the  various  parameters  needed  to  be  converted 
to  engineering  units.   This  was  done  to  make  it  easier  to 
program  the  calculations  necessary  to  compute  the  per- 
formance parameters  of  the  engine  and  to  put  the  data  in  a 
form  that  could  be  easily  understood  by  the  user. 

To  calibrate  the  digital-to-analog  converter,  data  was 
collected  for  a  number  of  points  within  the  range  that  was 
expected  to  be  measured  by  each  A/D  converter  channel.   A 
linear  curve  was  then  fit  to  the  data  and  the  curve  fit 
coefficients  entered  as  part  of  the  data  acquisition 
program.   The  raw  data  used  to  calibrate  each  A/D  converter 
channel  is  included  in  Appendix  B.   Plots  of  both  the  data 
and  the  linear  curve  fit  equations  are  included  in 
Appendix  C. 
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V.   DATA  ACQUISITION  PROGRAM 

The  data  acquisition  and  data  reduction  program,  listed 
in  Appendix  D,  was  written  entirely  in  FORTRAN  IV.   The 
program  was  written  on  a  Digital  VT-103  computer,  similar  to 
the  system  used  with  the  diesel  engine,  which  was  equipped 
with  an  8  inch  dual  disk  drive.   The  dual  disk  drive  allowed 
the  use  of  the  FORTRAN  compiler.   After  the  program  was 
written  and  compiled,  the  compiled  code  was  stored  on  a  tape 
cassette  using  one  of  the  two  internal  tape  drives  of  the 
VT-103*   The  tape  was  used  to  load  the  program  into  the 
computer  dedicated  to  data  acquisition  on  the  diesel  engine. 

There  were  two  sets  of  requirements  that  the  data 
acquisition  and  data  reduction  program  were  required  to 
meet.   One  set  of  requirements  determined  what  functions  the 
program  was  to  perform.   The  second  set  of  requirements  had 
to  do  with  making  the  program  "user  friendly". 

A.   PROGRAMMING  REQUIREMENTS 

The  basic  functions  of  the  program  were  to  read  a  number 
of  engine  parameters,  reduce  them  to  a  form  that  could  be 
displayed  to  the  user  and  use  the  recorded  parameters  to 
calculate  additional  quantities  used  in  performance 
evaluations  of  the  engine.   Additionally,  the  program  had  to 
provide  a  means  for  the  user  to  change  the  operating  point 
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of  the  engine  and  provide  a  means  of  entering  any  data 
necessary  for  the  calculations  that  could  not  be  measured 
directly  by  the  computer.   Finally,  the  program  had  to 
display  both  the  recorded  and  calculated  variables  and  print 
a  copy  of  the  display  if  required  by  the  user. 

In  order  to  meet  all  of  the  functional  requirements 
the  program  needed  to  be  written  to  run  interactively. 
Since  the  program  was  to  be  used  in  a  classroom  laboratory 
environment  it  needed  to  be  made  as  user  friendly  as 
possible.   It  would  be  used  by  students  who  would  have  vir- 
tually no  time  available  to  learn  how  the  program  operated. 
As  a  result,  it  had  to  be  almost  self-explanatory  when 
installed  and  running.   There  was  one  additional  require- 
ment, not  related  to  either  of  the  other  two  broad  sets  of 
requirements  that  needed  to  be  met  by  the  programmer. 

The  Digital  VT-103  computer  when  configured  as  the  one 
in  this  experiment  was,  has  an  available  memory  of  about  32K 
of  16-bit  words.   The  only  storage  devices  for  programs  were 
two  internal  cassette  tape  drives.   The  entire  program  had 
to  reside  in  the  computer  memory  that  was  available  after 
the  operating  system  was  installed  and  memory  partitioned 
for  variable  arrays  necessary  to  the  program.   Since  the 
drives  operate  very  slowly,  it  was  not  practical  to  have  a 
program  so  long  that  any  subroutine  would  need  to  be  called 
from  a  tape  while  the  program  was  running. 
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B.   PROGRAM  STRUCTURE 

The  data  acquisition  and  data  reduction  program  was 
written  as  a  set  of  subroutines,  each  designed  to  handle  a 
series  of  related  functional  requirements.   A  short  main 
program  was  written  to  allow  the  user  to  call  each 
subroutine  as  required  and  to  allow  the  user  to  exit  the 
program  if  desired. 

There  were  three  main  subroutines.   One  subroutine 
provided  a  means  of  inputting  the  values  for  four  variables, 
not  directly  measurable  by  the  computer,  these  were  needed 
in  the  data  reduction  section  of  the  program.   A  second 
subroutine  allowed  the  user  to  change  the  RPM  and  load  set 
points  of  the  engine.   The  third  subroutine  performed  the 
data  acquisition  and  data  reduction  functions  and  displayed 
and  printed  the  results. 

A  very  short  fourth  subroutine  was  added  after  the 
program  had  been  used  for  the  first  time  by  a  class  of 
students.   The  only  function  of  this  last  subroutine  was  to 
read  the  engine  RPM  and  torque  and  display  the  values  on  the 
terminal.   This  subroutine  was  provided  as  a  convenience  for 
the  operator. 

There  was  a  separate  menu  associated  with  the  main 
program  and  each  subroutine  that  listed  all  of  the  options 
available  to  the  user  while  in  each  subroutine.   Commands 
were-  issued  to  the  computer  by  typing  the  number  of  the 
option  desired. 
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C.   MAIN  PROGRAM 

When  the  computer  is  first  powered  up  there  is  no 
operating  system  that  automatically  boots  other  than  a 
simple  machine  language  monitor.   The  RT-11SJ  operating 
system  is  booted  from  a  tape  installed  in  one  of  the  tape 
drives.   After  installation  of  the  operating  system,  the 
data  acquisition/data  reduction  program  is  loaded  and  run 
from  the  other  tape  drive  by  typing  the  command  'RUN 
DIESEL' . 

The  main  program,  which  is  entered  when  the  program  is 
first  run,  contains  only  a  few  lines  of  FORTRAN  code.   Its 
only  function  is  to  provide  a  starting  point  for  user  inter- 
action with  the  program.   It  provides  a  means  of  calling  the 
subroutines  and  a  means  of  exiting  the  program. 

The  first  menu,  Figure  5,  to  appear  on  the  screen  is  the 
menu  associated  with  the  main  program.   The  program  uses  the 
prompt  '>'  anytime  it  is  ready  to  receive  a  command  through 
the  keyboard.   Commands  are  issued  by  simply  typing  the 
number  of  the  option  desired  followed  by  <return>.   Typing  1 
<return>  will  call  the  subroutine  that  samples  data  and 
displays  the  results  (subroutine  RUN).   Option  2  calls  the 
subroutine  that  allows  input  of  data  through  the  keyboard 
(subroutine  DATIN).   Typing  3  allows  the  user  to  change 
engine  RPM  or  load  (subroutine  ENGSET).   Typing  4  ends 
execution  of  the  program. 
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D.  SUBROUTINE  DATIN 

There  are  four  variables  used  in  the  data  reduction 
portion  of  the  program  that  are  not  measured  directly  by  the 
computer.   These  variables  must  be  entered  by  hand  each  time 
the  program  is  run.   There  are  no  default  values  for  these 
variables  other  than  0.   The  four  variables  are  atmospheric 
pressure,  specific  gravity  of  the  fuel,  the  temperature  of 
the  fuel  sample  and  the  diameter  of  the  air  inlet  nozzle. 
The  air  inlet  can  be  through  one  of  two  different  size 
nozzles.   When  running  the  program  this  subroutine  is 
usually  the  first  subroutine  called  since  subroutine  RUN 
cannot  be  called  until  this  data  has  been  input. 

Whenever  this  subroutine  is  entered,  the  menu  shown  in 
Figure  6  appears  on  the  computer  terminal.   Typing  1 
following  the  prompt  will  cause  the  computer  to  display, 
'Enter  new  value  for  atmospheric  pressure.'.   The  user  types 
in  the  new  value  and  the  computer  will  rewrite  the  original 
DATIN  menu  with  the  new  value  for  atmospheric  pressure 
displayed.   Typing  2,  3  or  4  will  result  in  a  similar 
response  by  the  computer.   Typing  5  will  exit  the 
subroutine . 

E.  SUBROUTINE  ENGSET 

The  ENGSET  subroutine  has  only  one  functiton,  it  allows 
the  user  to  choose  the  operating  point  of  the  engine.   To  do 
this,  the  program  must  accept  new  set  points,  convert  the 
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set  points  into  values  that  can  be  used  by  the  control 
computer  and  send  the  converted  values  to  the  control 
computer . 

When  entering  this  subroutine  the  following  menu  appears 
on  the  screen  (Figure  7).   Typing  1  or  2  allows  the  user  to 
change  the  current  values.   After  typing  one  of  the  first 
two  commands  the  computer  will  display,  'Input  an  integer 
value  for  desired  RPM  (or  torque).'.   The  user  can  then  type 
in  the  new  operating  point.   After  entering  the  new  point 
the  computer  will  display  the  first  menu  again  with  the  new 
value  for  either  RPM  or  torque  displayed.   The  computer  also 
takes  the  new  value  and  sends  it  to  the  control  computer. 

Before  sending  the  new  set  point  the  computer  must 
convert  it  to  a  value  that  can  be  passed  over  the  simple 
parallel  line  interface  connecting  the  two  computers.   The 
control  computer  uses  an  8-bit  microprocessor  while  the  data 
acquisition  computer  uses  a  16-bit  microprocessor.   A 
control  computer  'word'  is  eight  binary  bits  long  while  a 
data  acquisition  'word'  is  sixteen  binary  bits  long.   As 
mentioned  in  Chapter  3,  the  two  computers  are  interfaced  by 
connecting  the  eight  high  order  bits  from  the  VT-103 
parallel  port  to  one  of  the  eight  bit  parallel  ports  of  the 
control  computer  and  the  eight  low  order  bits  to  a  second 
parallel  port.   One  of  the  control  computer  parallel  ports 
receives  desired  torque  and  one  receives  desired  RPM,  but 
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the  16-bit  parallel  port  of  the  VT-103  has  to  send  both 
desired  torque  and  desired  RPM  at  the  same  time. 

An  8-bit  word  can  be  any  decimal  integer  between  0  and 
255.   A  16-bit  word  can  be  any  decimal  integer  between  0  and 
65536.   If  the  highest  bit  is  used  to  indicate  sign,  a 
16-bit  word  can  be  any  integer  between  -32768  and  +32768. 
The  VT-103  operating  system  uses  signed  numbers  so  the 
second  case  applied  in  this  instance.   Desired  torque  and 
RPM  are  translated  into  integers  between  0  and  32768  before 
being  sent  to  the  control  computer.   The  program  does  this 
by  first  converting  the  desired  RPM  to  an  integer  between  0 
and  255  using  equation  5.1. 


R  =  0.0755  *  (desired  RPM) 


(eqn .  5.1) 


Desired  torque  is  converted  to  an  integer  between  0  and  127 
using  equation  5.2. 


L  =  1.85  *  (desired  torque) 


(eqn.  5.2) 


The  integer  value  of  desired  torque  (L)  is  then  multiplied 
by  256.   This  has  the  same  effect  as  shifting  the  binary 
value  of  the  integer  to  the  eight  high  bits  of  a  16-bit 
word.   The  two  converted  values  are  then  added  together  and 
this  value  is  written  to  the  parallel  port  where  it  is 
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interpreted  as  two  separate  numbers  by  the  control  computer. 
The  integer  that  is  finally  written  to  the  parallel  port  is 
latched  by  both  computers.   This  means  that  the  number  is 
not  erased  until  another  number  is  written  to  the  port. 

Every  time  a  new  value  is  entered  by  the  user  for 
desired  RPM  or  desired  torque  the  program  rewrites  the 
number  in  to  the  parallel  port. 

Option  4  of  the  subroutine  command  options  allows  the 
user  to  exit  the  RPM  and  load  setting  subroutine.   Option  3 
calls  subroutine  SPEED.   This  subroutine  displays  actual 
RPM  and  actual  torque  and  then  returns  to  the  calling 
routine . 

F.   SUBROUTINE  RUN 

Subroutine  RUN  performs  the  most  functions  of  all  the 
routines  and,  as  a  result,  is  the  longest  and  most  complex 
portion  of  the  entire  program.   The  primary  functions  of  the 
RUN  subroutine  are  to  read  data  from  the  analog-to-digital 
(A/D)  converter  and  display  the  results.   This  subroutine 
performs  all  of  the  calculations  needed  to  convert  the 
output  from  the  A/D  converter  to  engineering  units  and  to 
compute  the  parameters  necessary  for  performance  evaluations 
of  the  engine.   Both  the  measured  parameters,  converted  to 
engineering  units,  and  the  computed  parameters  are  displayed 
on  the  terminal  and  can  be  printed  by  selecting  one  of  the 
command  options.   Other  command  options  allow  the  user  to 
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display  the  raw  output  from  the  ADC  with  the  standard 
deviation  computed  for  each  channel  or  to  call  the  ENGSET  of 
SPEED  subroutines. 

The  subroutine  RUN  menu  is  shown  in  Figure  8.   There  are 
seven  commands  that  can  be  given  while  this  subroutine  is 
executing . 

Command  option  1  will  cause  the  program  to  read  the 
first  12  channels  of  the  A/D  converter  and  performs  the 
calculations  necessary  to  convert  the  output  to  engineering 
units.   It  will  also  compute  the  performance  parameters; 
brake,  friction  and  indicated  horsepower;  mechanical  and 
thermal  efficiency;  mean  effective  pressure,  specific  fuel 
consumption,  scavenge  ratio,  air  mass  flow  rate,  and  air-to- 
fuel  ratio.   Appendix  E  explains  all  of  the  equations  used 
in  these  calculations.   After  performing  the  data 
acquisition  and  data  reduction,  the  program  will  display  the 
results  on  the  VT-103  monitor.   Figure  9  is  a  sample  of  the 
line-printer  output.   The  terminal  display  is  the  same 
except  that  it  includes  a  brief  form  of  the  command  menu 
displayed  below  the  output.   A  printed  copy  of  the  results 
can  be  obtained  by  issuing  command  2. 

Command  5  performs  the  same  functions  as  command  1  (data 
acquisition  and  reduction)  except  that  the  terminal  display 
is  changed.   The  screen  displays  the  mean  and  standard 
deviation  of  the  raw  data  from  each  channel  of  the  ADC.   The 
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command  1  display  can  be  obtained  by  issuing  command  3.   It 
is  not  possible  to  obtain  a  printed  copy  of  the  mean  and 
standard  deviations.   Command  5  was  used  to  calibrate  the 
transducer  output  to  the  A/D  converter. 

Command  4  calls  subroutine  ENGSET.   This  command  allows 
the  user  the  option  of  changing  the  operating  point  of  the 
engine  without  having  to  return  to  the  main  program  first. 
Command  7  calls  subroutine  SPEED  which  displays  the  actual 
RPM  and  torque.   Command  6  exits  the  subroutine. 
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VI.   CONTROL  THEORY 

Prior  to  designing  the  control  software  a  mathematical 
model  of  the  engine  and  dynamometer  transfer  functions  had 
to  be  developed  using  control  theory  and  experimental  data. 
Once  the  transfer  functions  were  determined,  appropriate 
control  actions  could  be  analytically  tested. 

A.   DETERMINATION  OF  TRANSFER  FUNCTIONS 

A  block  diagram  of  the  engine  and  dynamometer  is  shown 
in  Figure  10  with  the  measureable  parameters  being  the 
voltage  applied  to  the  throttle  (V) ,  throttle  position  (X) , 
engine  speed  (0),  voltage  to  the  dynamometer  armature 
current  controller  (I)  and  dynamometer  torque  (T).   Both  the 
engine  and  dynamometer  were  modeled  as  first  order  systems. 
There  were  two  additional  characteristics  of  the  engine  that 
would  conceivably  need  to  be  accounted  for  in  the  design  of 
the  control  system;  the  time  delay  (transportation  lag)  of 
the  engine  and  the  effect  of  the  throttle  position  control 
loop.   Because  of  this,  the  throttle  was  also  modeled  as  a 
first-order  system  and  an  exponential  term,  to  describe  the 
transportation  lag,  was  introduced  into  the  transfer 
function  for  the  engine.   Both  the  transportation  lag  and 
the  effect  of  the  transfer  function  of  the  throttle  were 
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eventually  determined  to  have  a  negligible  effect  on  the 
design  of  the  control  system. 

An  open-loop  first  order  transfer  function  has  the  form 
given  in  equation  6. 1 . 


9(3) 

ETsT 


K 

st  +  1 


(eqn .  6.1) 


where  K  is  the  steady  state  gain  constant  and  t  is  the  time 
constant  for  the  system.   Both  K  and  t  can  be  experimentally 
determined  for  most  systems.   To  determine  K  for  the  engine, 
the  dynamometer  load  was  set  to  zero  and  the  voltage  to  the 
throttle  was  recorded  for  various  engine  speeds.   During 
steady  state  operation,  K  =  e/V.   Figure  11  lists  the 
recorded  data  that  was  used  to  compute  K.   The  average 
magnitude  of  K,  with  units  of  ( 1 /volt-second) ,  for  the  14 
data  points  is  71-3-   This  value  was  used  in  the  initial 
analysis  of  the  control  problem. 

While  a  value  of  K  in  units  of  (1/volt-sec)  was 
convenient  for  the  initial  analysis,  once  the  computer  was 
installed  it  was  easier  to  compute  a  value  of  K  based  on 
the  decimal  number  output  to  the  digital-to-analog  converter 
and  use  the  new  value  of  K  in  further  analysis.   The  value 
of  the  open  loop  gain  constant  computed  by  dividing  engine 
RPM  by  D/A  value  is  13.2  RPM.   This  second  value  is  what 
appears  in  subsequent  calculations.   A  value  for  the  K  of 
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the  dynamometer  was  computed  in  the  same  manner.   The 
computed  gain  constant  for  the  dynamometer  is  equal  to 
0.48  lbf-ft. 

For  a  first  order  system,  "t"  is  equal  to  0.063  times 
the  settling  time  (t  )  of  the  system  after  a  step  input. 
For  the  diesel  engine  t.,  was  experimentally  calculated  by 
recording  the  response  of  the  engine  to  a  step  input. 
Figure  12  shows  a  representative  transient  response  curve  of 
the  engine  to  a  throttle  voltage  change  that  corresponds  to 
a  change  of  roughly  200  RPM.   A  2-pen  strip  chart  recorder 
was  used  to  record  the  transducer  outputs  for  throttle 
position  and  engine  speed.   The  transient  response  of  the 
engine  was  recorded  after  a  step  change  was  made  to  the 
throttle  voltage.   No  dynamometer  load  was  applied  to  the 
engine  for  these  measurements.   The  settling  time,  while 
found  to  vary  somewhat  depending  on  the  initial  engine 
speed,  was  determined  to  be  approximately  120  seconds.   This 
gave  a  time  constant  for  the  engine  of  approximately  75 
seconds . 

The  time  constant  (T)  associated  with  the  transportation 
lag,  while  difficult  to  determine  experimentally  without  the 
benefit  of  more  sophisticated  electronic  equipment  can  be 
approximately  determined  analytically.   Neto  and  Thaler 
[Ref.  6]  showed  that  T  could  be  approximated  for  a  diesel 
engine  by  equation  6.2. 
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T  =  40/RPM 


(eqn.  6.2) 


The  operating  range  of  the  General  Motors  diesel  was  1500  to 
2500  RPM.   These  speeds  corresponded  to  time  constants  for 
the  exponential  term  between  0.016  and  0.027  seconds.   As  a 
result,  it  could  be  assumed  that  the  effects  of  transpor- 
tation lag  only  become  significant  at  frequencies  higher 
than  the  operating  range  of  the  engine  and  that  the  effect 
of  transportation  lag  can  be  neglected. 

The  time  constant  for  the  first  order  model  of  the 
throttle  positioning  loop  was  determined  using  the  same 
technique  as  used  to  determine  the  time  constant  for  the 
diesel  engine.   Settling  time  for  the  throttle  positioning 
loop  was  less  than  0.5  seconds.   This  gave  a  time  constant 
for  the  throttle  (tp)  of  less  than  0.315  seconds.   As  in  the 
case  of  transportation  lag,  the  response  time  of  the 
throttle  positioning  loop  can  be  ignored  since  it  becomes 
significant  only  at  frequencies  higher  than  the  operating 
frequencies  of  the  engine. 

The  steady-state  gain  constant  for  the  throttle  cannot 
be  ignored,  but,  because  data  was  recorded  for  throttle 
voltage  vs.  RPM  no  additional  calculations  are  required. 
The  13-2  RPM  gain  constant  is  in  actuality  the  throttle 
gain  constant  multiplied  by  the  engine  gain  constant.   The 
final  diesel  engine  transfer  function  for  speed  of  the 
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engine  as  a  function  of  the  input  to  the  throttle  is  given 
in  equation  6.3. 


RPM(s)  13.2e~,Q2s  ,  c    ,, 

E^TT     =      (75s   +    1X.315S   +    1)  (eqn'    6'3) 


The  Bodie  diagrams  for  the  open  loop  transfer  function  of 
the  diesel  as  determined  experimentally  is  shown  in  Figure 
12  and  Figure  13  curve  1. 

The  same  analysis  was  used  to  determine  the  transfer 
function  of  the  dynamometer.   For  the  dynamometer  the 
armature  current  could  be  varied  to  control  the  torque 
applied  to  the  engine.   The  magnitude  of  the  load  was  a 
function  of  the  voltage  applied  to  the  armature  current 
controller.   The  open  loop  gain  constant  was  experimentally 
determined  at  a  constant  engine  speed  of  1500  RPM. 

The  steady  state  dynamometer  gain  can  be  shown  to  be  a 
function  of  engine  RPM  but  for  the  purpose  of  this  analysis 
was  assumed  to  be  a  constant,  0.48  Ibf-ft.   The  dynamometer 
had  a  calculated  time  constant  that  was  very  much  smaller 
than  the  diesel  engine.   A  plot  of  the  transient  response  of 
the  open  loop  transfer  function  was  obtained  in  a  manner 
similar  to  the  method  used  to  determine  the  time  constant 
for  the  engine.   A  2-pen  strip  chart  recorder  was  connected 
to  the  voltage  supply  to  the  armature  current  controller  and 
the  dynamometer  load  sensor.   The  dynamometer  had  a  settling 
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time  of  approximately  1.0  second  which  gave  a  time  constant 
(tO  of  0.63  seconds.   The  open  loop  transfer  function  for 
the  dynamometer  as  determined  experimentally  is  given  in 
equation  6.4. 


T(s)  0.48  ,    nn      ,    us 

E^riT  =   (.63s  +  1)  (eqn*  6'4) 

The  Bodie  plots  of  the  dynamometer  transfer  function  are 
shown  in  Figures  14  and  15  curve  1. 

B.   DESIGN  OF  CONTROL  ACTION 

In  designing  the  control  action  for  both  the  engine  and 
dynamometer  the  analysis  was  made  assuming  that  the  two 
systems  were  independent.   A  proportional-plus-integral- 
plus-differential  control  action  was  chosen  for  the 
controller.   The  use  of  a  microprocessor  allowed  the  use  of 
any  control  action  since  the  action  depended  only  on  the 
algorithm  programmed  into  the  computer  and  not  on  the  system 
hardware.   The  PID  controller  offered  several  character- 
istics that  were  desirable  in  a  control  system. 

A  proportional  controller  has  a  function  in  the  time 
domain  of  as  shown  in  equation  6.5. 

e(t)  =  Ae(t)  (eqn.  6.5) 

In  the  Laplace  domain  the  function  was  transformed  into 
equation  6.6. 
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|{f}   =   A  (eqn.  6.6) 


The  A  in  each  function  is  an  adjustable  gain  constant  and 
allowed  the  programmer  a  method  of  adjusting  the  sensitivity 
of  the  control  system. 

The  equation  for  an  integral  control  action  in  the  time 
domain  is  given  in  equation  6.7  for  the  time  domain  and  in 
equation  6.8  for  the  Laplace  domain. 


C 

\      e(t)dt 


e(t)   =   B   J   e(t)dt  (eqn.  6.7) 


e(s)     B  ,     ,  QS 

EUT     =  i"  (eqn*  6*8) 


An  integral  control  action  has  the  characteristic  of 
providing  essentially  zero  steady  state  error  to  a  control 
system.   The  integral  of  a  small  error  between  the  actual 
operating  point  and  the  desired  operating  point  will 
eventually  become  large  enough  to  cause  the  system  to  return 
to  the  desired  operating  point.   The  "B"  in  the  transfer 
function  is  a  programmable  gain  associated  with  the  integral 
control  action . 

A  differential  control  action  will  respond  only  to  a 
change  in  error  and  as  a  result  can  only  be  used  in 
conjunction  with  a  proportional  or  integral  control  action. 
It  has  the  unique  characteristic  of  only  effecting  a 
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response  if  the  error  in  the  controlled  parameter  is 
changing.   This  characteristic  makes  it  impossible  to  use 
only  a  differential  control  action  in  a  control  system  but 
when  used  in  conjunction  with  one  or  both  of  the  other 
control  actions  the  controller  can  anticipate  the  approach 
of  the  system  to  the  desired  control  point  and  slow  the 
action  of  the  controller.   The  function  for  a  differential 
control  action  has  the  function  in  the  time  domain  as  shown 
in  equation  6.9  and  a  function  in  the  Laplace  domain  as 
shown  in  equation  6.10. 


e(t)   =   C  4|i£!  (eqn.  6.9) 


|i|l   =   Cs  (eqn.  6.10) 


The  constant  C  is  a  variable  gain  constant  associated  with 
the  differential  control  action. 

One  additional  approximation  is  required  when  a  digital 
computer  is  installed  in  the  feedback  control  loop.   A 
digital  computer  does  not  compute  the  error,  error  integral 
and  error  differential  continuously.   The  computer  samples 
once  during  a  discrete  time  period.   In  this  application, 
the  computer  samples  once  every  0.54  seconds;  the  length  of 
time  it  takes  to  complete  one  loop  in  the  control  program. 
The  continuous  functions  for  the  integral  and  differential 
control  actions  are  approximated  by  equation  6.11  and  6.12. 
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t 

e(t)   =   BT_       e„  (eqn.  6.11) 

n  =  0 

e(t)   =  J-   (en  -  en-1)  (eqn.  6.12) 

s 

Combining  the  three  control  actions  gives  an  's'  domain 
function  as  given  in  equation  6.13. 


e(s)  BTS     C 

EXIT  =  A  +  "S"  +  —  S  (eqn'  6'13) 

A  small  amount  of  algebraic  manipulation  puts  the  function 
in  Bodie  form  (equation  6.14). 


T  B  (  -£-x  S2  +  —  S  +  1)       (eqn.  6.14) 
e(s)  BTS^      BTS 

ElsT   =  S 


The  steady  state  gain  constant  for  the  function  in  Bodie 
form  now  becomes  B*TS.   The  natural  frequency  associated 
with  the  controller  ( u>)  is  equal  T  *sqrt(B/C)  and  the 
damping  ratio  for  the  controller  (5)  is  equal  to, 


A/sqrt(4*B*C)  (eqn.  6.15) 


The  T   in  each  equation  represents  the  sample  time  period; 
the  time  the  control  computer  takes  to  make  one  iteration  of 
the  control  loop.   For  this  system  T   =  0.54  seconds. 

S 
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The  Bodie  diagram  for  the  diesel  engine  PID  controller 
are  shown  in  Figure  12  and  Figure  13  curve  2.   The  Bodie 
plots  for  the  dynamometer  are  given  as  curve  2  of  Figures  14 
and  15. 

The  initial  design  of  the  controller  necessitated 
determining  values  for  the  three  constants  associated  with  a 
PID  controller  (A,  B  and  C) .   These  constants  were 
originally  determined  analytically  but  their  final  values 
were  obtained  experimentally  by  monitoring  the  performance 
of  the  system  with  the  control  computer  installed. 

The  open  loop  Bodie  plots  for  the  compensated  diesel 
engine  are  shown  in  Figures  12  and  13  curve  3.   These  are 
the  diagrams  of  the  system  in  its  final  form.   The  objective 
in  designing  the  system  was  to  provide  an  open  loop  gain 
margin  of  at  least  5  dB  and  a  phase  margin  of  approximately 
45  degrees.   In  both  cases  the  margins  were  exceeded.   The 
final  values  for  the  three  constants  for  the  diesel  engine 
speed  control  were  A=1,  B=0.0125,  and  C=5.   For  the  dyna- 
mometer load  controller  the  values  were  A=0.04,  B=0.004,  and 
C=0.2.   As  stated  previously,  Figures  12  and  13  show  the 
Bodie  plots  of  the  uncompensated  and  compensated  diesel 
engine  and  the  diesel  engine  controller  transfer  functions. 
A  block  diagram  for  the  controlled  diesel  engine  and 
dynamometer  showing  all  of  the  transfer  functions  is  shown 
in  Figure  16.   The  Bodie  plots  for  the  compensated 
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dynamometer  transfer  function  are  given  in  Figures  14  and  15 
curve  3  along  with  the  uncompensated  dynamometer  and  the 
dynamometer  controller  transfer  functions. 
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VII.   CONTROL  SOFTWARE 

The  control  software  that  is  installed  in  the  SYS-1 
computer  is  written  entirely  in  interpreted  NSC  Tiny  BASIC. 
The  control  algorithm  is  very  easy  to  program  in  BASIC  and 
requires  only  a  few  lines  of  code.   In  the  final  program 
limits  are  placed  on  the  numerical  value  of  certain  program 
variables  to  prevent  register  overflow  and  to  lessen  the 
impact  of  system  noise  on  the  response  of  the  system 

The  basic  coding  of  the  control  algorithm  requires  only 
7  lines  of  programming. 

10  E  =  R-M 
20  Isl+E 

30  D=E-E1 

40  E1=E 

50  V=A*E+B*I+C*D 

60  N=N+V 

70  GOTO  10 

In  the  above  example,  the  error  (E)  is  determined  at  line  10 
by  subtracting  the  measured  value  (M)  for  a  controlled 
parameter  from  the  requested  value  (R).   In  line  20  a 
numerical  integration  is  performed  by  adding  the  current 
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error  to  the  error  total  from  all  previous  iterations.   Line 
30  performs  a  numerical  differentiation  by  subtracting  the 
error  from  the  previous  iteration  from  the  current  error. 
The  numerical  values  for  the  error,  error  integral  and  error 
differential  are  multiplied  by  their  respective  gain 
constants  (A,  B  and  C)  and  added  together  to  give  a 
numerical  value  representing  the  change  (V)  in  the  output  to 
the  actuator  (N).   In  line  60  the  change  is  added  to  the 
existing  value  and  line  70  returns  the  program  to  line  10  to 
begin  another  iteration.   The  actual  control  program  is  not 
quite  so  simple.   All  further  comments  on  the  control 
program  refer  to  the  listing  of  the  actual  control  program 
found  in  Figure  17. 

The  explanation  of  the  program  installed  in  the  SYS-1 
begins  at  line  100.   Actual  engine  speed  and  load  are  input 
through  channels  5  and  8,  respectively,  of  the  analog-to- 
digital  converter.   Channel  5  is  located  at  address  #9804 
and  channel  8  is  located  at  address  #9808.   NSC  Tiny  BASIC 
supports  most  of  the  commands  found  in  Dartmouth  BASIC  but 
many  require  a  shortened  code.   A  listing  of  the  format  for 
Tiny  BASIC  commands  can  be  found  in  the  SYS-1  User's  Guide 
[Ref.  3].   In  addition  to  only  using  Tiny  BASIC  as  a  higher 
level  language,  the  SYS-1  can  only  perform  integer 
arithmetic . 
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The  @  command  in  Tiny  BASIC  performs  the  same  functions 
as  PEEK  and  POKE  in  standard  BASIC.   The  '#'  sign  indicates 
a  number  is  a  hexidecimal  number.   The  explanation  of  the 
code  begins  at  line  100  which  reads  the  actual  values  for 
engine  RPM  and  engine  load.   The  analog-to-digital  converter 
must  be  written  to  in  order  to  begin  a  conversion,  so  0  is 
written  to  addresses  #9804  and  #9808  prior  to  reading  from 
addresses  #9804  and  #9808.   The  ADC  requires  8  clock  cycles 
to  perform  one  conversion.   This  is  less  time  than  it  takes 
the  microprocessor  to  interpret  the  next  command  so  there  is 
no  danger  of  attempting  to  read  a  channel  prior  to  the 
completion  of  a  conversion . 

Requested  values  for  engine  RPM  and  torque  are  read  at 
line  110.   Addresses  #0A00  and  #0A01  are  the  locations  for 
Port  A  and  Port  B  of  the  SYS-1  computer's  installed 
peripheral  interface  adapter  (PIO).   These  two  ports  are 
used  to  interface  the  SYS-1  to  the  parallel  line  interface 
of  the  VT-103  data  aquisition  computer.   Due  to  the 
configuration  of  the  VT-103  operating  system,  the  largest 
integer  number  that  can  be  transferred  through  Port  B  is 
127.   Because  of  this,  the  value  received  at  register  #0A00 
is  multiplied  by  3  to  increase  its  value  to  the  same  range 
as  that  read  at  register  #9808. 

If  the  value  for  requested  RPM  is  less  than  10  (10 
corresponds  to  approximately  100  RPM) ,  the  program  jumps  to 
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line  10  for  execution.   Lines  10,  20  and  30  form  a  secondary 
loop  outside  the  main  program  that  stops  execution  of  the 
control  algorithm  until  a  value  greater  than  10  is  input 
through  the  PI0.   This  secondary  loop  smooths  the  transition 
from  manual  to  computer  control  of  the  engine  and  dyna- 
momter.   Prior  to  returning  to  the  main  program  the  error 
integrals  are  re-initialized  and  the  current  engine  RPM  is 
read  and  the  value  for  engine  RPM  is  written  to  the  register 
controlling  throttle  position. 

Lines  120  through  170  determine  the  control  action  of 
the  computer  and  output  to  the  actuators.   There  are  10 
IF-THEN  statements  within  this  portion  of  the  program  which 
serve  to  limit  the  magnitude  of  the  absolute  value  of 
certain  variables.   These  statements  prevent  some  variables 
from  overriding  others  and  producing  adverse  effects  in  the 
control  action  and  also  prevent  overflow  errors  from 
occurring  in  the  registers  associated  with  the  digital-to- 
analog  (D/A)  converters. 

Limits  are  placed  on  the  error  integrals  and  output  for 
both  the  RPM  and  torque  actuators.   In  addition  a  limit  is 
placed  on  the  magnitude  of  the  change  that  can  be  output  to 
the  throttle  position  between  any  two  successive  iterations 
of  the  engine  speed  control  algorithm. 

The  limit  on  V,  the  change  in  throttle  position, 
protects  the  system  from  spurious  conversions  by  the  ADC  due 
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to  electrical  noise.   During  system  operation,  it  was 
noticed  that  occasionally  inconsistent  readings  would  be 
received  by  the  computer  from  the  ADC.  The  engine  would  be 
running  at  a  steady  RPM  and  suddenly,  during  one  iteration, 
the  computer  would  read  an  error  of  several  hundred  RPM. 
The  error  invariably  returned  to  a  value  consistent  with 
previous  iterations  on  the  following  iteration  but  the 
transient  behavior  of  the  system  due  to  the  spurious  reading 
could  be  quite  unsatisfactory,  requiring  a  minute  or  more  to 
return  to  steady  state,  if  the  magnitude  of  V  were  not 
limited . 

The  value  of  the  integral  term  in  the  control  equation 
is  greatest  when  the  system  is  operating  close  to  its 
desired  control  point.   The  integral  term  always  returns  the 
system  to  the  precise  control  point  requested  within  the 
accuracy  the  controller  can  measure.   The  disadvantage  of 
the  integral  term  occurs  when  the  system  is  slow  to  respond 
to  a  change  to  the  control  point.   As  an  example,  the  diesel 
engine  will  react  slowly  to  a  speed  change  if  it  is 
operating  under  a  heavy  load.   In  this  case  if  the  error 
integral  is  not  limited  it  can  become  very  large  and  the 
system  will  steady  out  at  a  point  above  or  below  the  control 
point  until  the  magnitude  of  the  error  integral  is  reduced 
to  near  zero.   Limiting  the  magnitude  of  the  error  integrals 
prevents  this  type  of  undesirable  system  response. 
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The  last  set  of  variables  artifically  limited  by  the 
program  are  the  variables  output  to  the  actuators.   These 
variables  are  written  to  a  PIO  which  sends  the  digital 
equivalent  of  the  variable  to  the  D/A  converters  to  be 
converted  into  analog  voltages.   The  value  of  the  variables 
written  to  the  D/A  converters  is  limited  to  integers  between 
0  and  255.   The  8-bit  DACs  used  in  this  system  cannot 
convert  integers  outside  this  rangee. 

Line  180  prints  the  current  values  for  all  variables 
prior  to  the  program  returning  to  line  100  to  state  another 
iteration.   While  this  line  is  not  necessary  for  the 
operation  of  the  system  it  aids  in  designing  and  fine  tuning 
the  program  and  presents  a  visual  representation  of  the 
transient  behavior  of  the  system. 
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VIII.   RESULTS  AND  CONCLUSIONS 

The  integrated  data  acquisition  and  data  reduction 
system  and  engine  control  system  meet  all  of  the  system 
criteria  but  both  systems  could  benefit  from  further 
development.   After  manually  starting  the  engine,  speed  and 
load  control  could  be  passed  to  the  control  computer  and  the 
entire  range  of  experiments  run  in  class  could  be  performed 
using  the  data  acquisition  computer.   The  operator  was  able 
to  change  engine  RPM  and  load,  measure  data  and  output  data 
using  the  computer  keyboard.   Improvements  could  be  made  in 
two  areas  of  each  system  (control  system  and  data  acquisi- 
tion system) .   Improvements  in  resolution  and  transient 
response  could  be  made  to  the  control  system  and  improve- 
ments in  data  storage  and  in  the  calculation  of  one  of  the 
performance  parameters  of  the  engine  could  be  made  in  the 
data  acquisition  system. 

A.   CONTROL  SYSTEM 

The  diesel  engine  control  system  will  only  maintain 
engine  RPM  and  load  within  the  resolution  of  the  8-bit 
analog-to-digital  converter.   The  maximum  theoretical 
resolution  of  the  converter  in  RPM  is, 

+(1/255)  x  2500  =  +9.8  RPM  (eqn.  8.1) 
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This  resolution  could  be  obtained  if  the  RPM  sensor  output 
+5  volts  when  the  engine  was  operating  at  its  maximum  speed 
(2500  RPM).   The  actual  measured  resolution  is  +13.1  RPM. 
As  a  result  of  the  poor  resolution  of  the  engine  speed,  the 
controlled  engine  speed  will  fluctuate  around  the  desired 
RPM.   This  characteristic  of  the  control  system  appears 
aggravated  when  the  engine  is  operating  at  high  speeds. 

The  problem  could  be  alleviated  by  using  a  12-bit  A/D 
converter  in  place  of  the  8-bit  converter.   A  12-bit 
converter  would  have  a  theoretical  resolution  of, 


+(1/4096)  x  3000  =  +0.6  RPM 


(eqn.  8.2) 


This  resolution  would  allow  much  more  precise  control  of 
engine  speed. 

Using  the  8-bit  converter  to  control  engine  load  does 
not  create  the  same  problems  as  it  does  with  speed  control. 
The  theoretical  resolution  of  the  8-bit  A/D  converter  when 
measuring  torque  is  +0.2  lbf-ft.   The  actual  measured 
resolution  is  +0.5  lbf-ft.  which  is  well  within  the 
acceptable  error  in  torque  for  the  experiments. 

The  second  characteristic  of  the  control  system  that 
would  benefit  from  further  development  is  the  transient 
response  of  the  controlled  system.   The  engine  exhibits  an 
overshoot  of  approximately  10%  when  responding  to  a  step 
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change  to  desired  RPM.   This  overshoot  occurs  both  when  the 
engine  is  accelerating  or  decelerating.   The  transient 
response  characteristics  of  the  controlled  system  could  be 
improved  by  optimizing  the  values  of  the  three  constants 
associated  with  each  controller. 

A  more  accurate  mathematical  model  of  the  transfer 
functions  for  the  engine  and  dynamometer  would  aid  in 
optimizing  the  control  system.   In  the  transfer  functions, 
the  steady  state  gain  of  each  system  is  not  a  constant. 
The  gain  of  the  engine  is  a  function  of  torque  (K  =  k(T)) 
and  the  gain  of  the  dynamometer  is  a  function  of  RPM  (K  = 
k(RPM)).   A  refined  form  of  the  diesel  engine  transfer 
function  is  given  in  equation  8.3  and  a  refined  form  of  the 
dynamometer  transfer  function  is  given  in  equation  8.4 


RPM(s)  k(T)e~TS  ,      n  B    -. 

•= — -f — r—        =  7— -v  ;■  ■ — - — ; rr  V  eQn  •  O.j 

E.  (  S)  (t-jS    +    1  )  (  toS     +     1  ) 

T(s)         _  k(RPM)  ,nn  a    Uv 

l^TsT       =  (t23     ♦     D  (eqn*  8*4) 


In  the  two  equations,  k(T)  and  k(RPM)  can  be  experimentally 
determined . 

B.   DATA  ACQUISITION  SYSTEM 

The  biggest  problem  with  the  data  acquisition  and  data 
reduction  system  is  due  to  the  lack  of  data  storage 
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capacity.   The  VT-103  tape  drives  operated  so  slowly  that  it 
was  impractical  to  store  data  on  tape  when  running  experi- 
ments.  The  system  needed  a  disk  drive  for  faster  data 
storage  and  access.   The  slow  speed  of  the  tape  drives  also 
caused  undesirable  delays  when  running  the  program.   Using 
the  installed  tape  drives  required  ten  minutes  to  boot  the 
operating  system  and  load  and  run  the  data  acquisition 
program.   Anytime  the  user  required  a  printout  of  the  data, 
a  delay  of  2.5  minutes  was  experienced  while  the  program 
called  a  system  subroutine  from  the  system  library  stored  on 
tape.   A  disk  drive  installed  in  the  system  would  have 
alleviated  all  of  these  problems. 

The  second  area  for  future  development  in  the  data 
acquisition  system  is  in  the  computation  of  engine 
frictional  horsepower.   This  is  the  horsepower  the  engine 
expends  in  overcoming  internal  friction  and  in  powering 
auxiliary  equipment  such  as  the  oil  pump.   The  program 
calculated  friction  horsepower  by  using  a  curve  fit  to  the 
results  of  prior  experiments.   The  data  for  the  curve  fit 
was  obtained  by  rotating  the  engine  with  the  throttle  closed 
using  the  dynamometer  as  an  electric  motor  and  recording  the 
generated  torque.   While  this  curve  appears  reliable  under 
most  conditions  it  can  give  unreliable  results  if  engine 
characteristics,  such  as  operating  temperatures,  change.   A 
subroutine  needs  to  be  written  which  will  power  the  engine 
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with  the  dynamometer  and  record  the  torque  as  a  function  of 
engine  RPM.   The  program  could  then  fit  a  curve  to  the 
acquired  data  that  would  give  a  more  accurate  representation 
of  engine  friction  horsepower  for  any  given  series  of 
experiments . 

A  last  comment  is  needed  on  the  calibration  of  the 
analog-to-digital  converter  channels.   Because  of  time 
constraints,  accurate  calibrations  were  not  made  for  any  of 
the  A/D  channels.   Instead,  the  transducers  were  left 
installed  and  the  values  read  from  the  converter  channels 
were  compared  to  the  values  read  on  the  existing  gages.   A 
range  of  values  was  recorded  for  each  channel  by  operating 
the  engine  through  a  variety  of  speed  and  load  conditions. 
The  readings  recorded  from  the  existing  gages  could  not  be 
as  accurate  as  desired  because  the  recorded  engine 
parameters  tend  to  fluctuate  while  the  engine  is  operating 
at  any  given  operating  point.   A  more  accurate  calibration 
of  the  analog-to-digital  converters  is  desired. 
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Figure  1.   Machinery  Layout;  Engine,  Dynamometer, 
and  Transducers 
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DIESEL  ENGINE  DATA  AQUISITION  AMD  CONTROL  PROGRAM 


C u r r e n t  c  o mm a n d  optio n s  are" 

1  .  Data  acquisition  and  display  s  u  b  r  o u t i  n  e 

2  .  I  Tip  u  t  d  a  t  a  thro  u  s  h  K  a  v  board  . 

3 .  E n 3 i  n  e  R PM  and  I o a d  s e  1 1  i n 3  subro u t  i n e . 

4.  Exit  prosram. 


Enter  desired  option. 


Figure  5.   Data  Acquisition  Program  Main  Command  Menu 
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r 

ENGINE  SPEED  AND  LOAD  SET 

C  Lirr  en  t  values!    R  PM  -     0    TORQUE  =   Of  t .  1  b 
C  u r r e n  t  c  o  m  m  and  options  are" 

1  =  to  c h a n a e  ansine  R  P  N 

2.  to  chanse  e n a i n e  load 

3 .  to  retur n 

4.  to  display  RPM  and  load 

Enter  desired  option. 


Figure  7.   Subroutine  ENGSET  Command  Menu 
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Engine 

0 

V 

K 

RPM 

(rad/sec ) 

(volts) 

( 1 /volt-sec) 

510 

53.41 

0.755 

70.7 

1000 

104.72 

1.471 

71.2 

1500 

157.08 

2.260 

69.5. 

2010 

210.49 

2.960 

71.  1 

2300 

240.86 

3.  180 

75.8 

2490 

260.75 

2.670 

71.0 

2000 

209.44 

2.940 

71  .2 

1520 

159. 17 

2.230 

71.4 

1500 

157.08 

2.210 

71  .  1 

1230 

128.81 

1.873 

71.0 

1000 

104.72 

1.510 

69.4 

920 

96.34 

1.358 

70.9 

550 

57.59 

0.778 

74.0 

Average  value  of  K  =  71.4 
Determined  at  no  load  conditions. 


Figure  11.   Data  for  Determination  of  Engine  Gain  Constant 
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1  REM:  DIESEL  ENGINE  CONTROL  PROGRAM 

5  @#980B=#80 

10  PRINT  "CONTROL  PROGRAM" 

20  R=@OAOO 

25    @#9804  =  0:N  =  @//9804 

30    IF    (R<10)    GOTO    20 

40    §#9808=N:I=0:A=0 

100  g#9804=0:M@#9804: §#9807=0 :Xs§#9807 

110  R=§#0A00:T=3*@#0A01 

120  IF  (R<10)  GOTO  10 

130  E=R-M:Y=T-X:I=I+E:A=A+Y:D=E-F:G=Y-Z 

131  V=E+I/80+5*D 

131  IF  (V>50)  THEN  V  =  50 

132  IF  (V<-50)  THEN  V=-50 
140  J=(Y+A/20+5*G)/25 

145  IF  (I>200)  THEN  1=200 

146  IF  (K-200)  THEN  I  =  -200 

147  IF  (A>150)  THEN  A=150 

148  IF  (A<-150)  THEN  A=-150 
150  F=E:Z=Y:N=N+V:L=L+J 

160  IF  (N>255)  THEN  N=255 

161  IF  (N<0)  THEN  N=0 

162  IF  (L>255)  THEN  L  =  255 

163  IF  (L<0)  THEN  L=0 
170  @#9808=N:@#9809=L 

180  PRINT  E,I/80,5*D,V,N,"      " , Y, A/20 , 5*G , J , L 
190  GOTO  100 
200  END 


Figure  17.   Control  Program  for  SYS-1  Computer 
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APPENDIX  A 
EQUIPMENT  AND  INSTRUMENTATION 

Table  A. 1 


Engine:  GM  Model  1-53X3 

2  stroke  cycle 
Roots  3  lobe  blower 
30  BHP  at  3000  RPM 


Dynamometer:  Reliance  Super  T  DC  motor 

Frame  366-A-DZ,  30  HP 
Serial  number  U  302715-71 


RPM:  Anadex  Panel  Meter  Counter 

Model  CPM-603,  ser .  43911 
0  to  99999  RPM 
note  1 . 


Torque:  Wallace  and  Tiernan  Model  FA  145 

Serial  number  DD  08995 
Scale  0-70  and  70-125  lbf-ft 
note  2. 


Exhaust  Gas  Pressure:   Celesco  LCVR  Differential  Pressure 

Transducer  with  LCCD  Carrier 
Demodulator 
Range  0-25  cm  H20 
0-5  volts  DC  output 


Inlet  Air  Pressure:     Celesco  LCVR  Differential  Pressure 

Transducer  with  LCCD  Carrier 
Demodulator 
Range  0-25  cm  H20 
0-5  volts  DC  output 


Boost  Air  Pressure:     Celesco  LCVR  Differential  Pressure 

Transducer  with  LCCD  Carrier 
Demodulator 
Range  0-100  cm  H20 
0-5  volts  DC  output 
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Fuel  Mass  Flow  Rate: 


Anadex  Panel  Meter  Counter 
Model  CPM-603,  ser .  43913 
note  1 . 


Fuel  Temperature: 


Newport  Pyrometer 
Model  268TF1 ,05,CC2 
Serial  99180149-25 
Scale  0-999.99  deg  F 
Output  0.001  volts/deg  F 


Exhaust  Gas  Temperature:  Newport  Pyrometer 

Model  268TF1,05,CC2 
Serial  9108656-25 
Scale  0-999.99  deg  F 
Output  0.001  volts/deg  F 


Inlet  Air  Temperature: 


Newport  Pyrometer 
Model  268TF1,05,CC2 
Serial  9180154-25 
Scale  0-999.99  deg  F 
Output  0.001  volts/deg  F 


Water  Temperature  In: 


Newport  Pyrometer 
Model  268TF1 ,05,CC2 
Serial  9180150-25 
Scale  0-999-99  deg  F 
Output  0.001  volts/deg  F 


Water  Temperature  Out: 


Newport  Pyrometer 
Model  268TF1,05,CC2 
Serial  9180151-25 
Scale  0-999.00  deg  F 
Output  0.001  volts/deg  F 


Note  1 .  -  The  Anadex  meters  used  to  measure  RPM  and  fuel  mass 
flow  have  variable  frequency  output  signals.   Instrumentation 
amplifiers  were  constructed  locally  to  convert  the  frequency 
output  to  voltage  output  in  the  range  of  0-5  volts  DC. 
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Note  2.  -  The  Wallace  and  Tiernan  torque  meter  measures  torque 
by  sensing  the  air  pressure  developed  by  a  bellows  attached  to 
a  lever  arm  that  positions  the  dynamometer.   A  Celsco  LCVR 
differential  pressure  transducer  and  LCCD  carrier  demodulator 
were  inserted  in  parallel  with  the  torque  meter  to  sense  the 
bellows  air  pressure.   Output  from  the  LCVR  was  converted  to 
lbf-ft  torque  by  comparing  it  to  the  existing  torque. 
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APPENDIX  B 
DATA  FOR  CALIBRATION  OF  THE  ANALOG- TO- DIGITAL  CONVERTER 


Table  B.l.   Channels  1-5,  Temperature 


D  i  g  i  ta  1   n    c 

u  i        De9  F 
Value       3 


Resolution  =  2.5  deg  F 
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APPENDIX  C 
CALIBRATION  CURVE  PLOTS 
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APPENDIX  D 
DATA  ACQUISITION  PROGRAM  LISTING 
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APPENDIX  E 
DATA  REDUCTION  CALCULATIONS 


1.   Fuel  Mass  Flow  Rate  (MFUEL); 

-  fuel  mass  flow  rate  read  directly  from  mass  flow 
meter . 

MFUEL  =  MF  (channel  9  of  A/D  converter) 


2.  Brake  Horsepower  (BHP);  net  rate  of  work 

BHP  =   (2  *  PI  *  RPM  *  TORQ)/33000 

3.  Friction  Horsepower  (FHP); 

-  compute  friction  torque  (FTORQ)  using  curve  fit 
equation  developed  from  experimental  data. 

FTORQ  =  10.636  +  (RPM  *  1.219E-02) 

-  compute  friction  horsepower  (FHP) 

FHP  =  (2  *  PI  *  RPM  *  FTORQ)/33000 

4.  Indicated  horsepower  (IHP);  total  engine  work  rate 

IHP  =  BHP  +  FHP 

5.  Mechanical  Efficiency  of  Engine  (MECH);  net  work  rate 
divided  by  total  work  rate 

MECH  =  BHP/IHP 

6.  Mean  Effective  Pressure  (MEP); 

MEP  =  24  *  PI  *  T  /  VD 
VD  =  displacement  volume  =  53*0685  in 

7.  Specific  Fule  Consumption  (SFC);  rate  of  fule  usage 
divided  by  net  work  rate 

SFC  =  (3600  *  MFUED/BHP 
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8.  Thermal  Efficiency  of  Engine  (THERM); 

-  convert  fuel  sample  specific  gravity  (SPGR)  to 
specific  gravity  at  60  deg  F  (SG60) 

SG60  =  SPGR  +  (FST-60)/3600 

-  compute  lower  heating  value  of  fuel  (LHV) 

LHV  =  16380  +  (API) 
API  =  141.5/SG60  -131.5 

-  compute  thermal  efficiency 

THERM  =  (0.70694  *  BHP)/(MFUEL  *  LHV) 

9.  Air  Mass  Flow  Rate  (MAIR)  through  air  inlet  nozzle 

Y  =  (PA  -  PAIR)/PA 
PA  =  atmospheric  pressure  (stagnation  pressure) 
PA  -  PAIR  =  pressure  differential  across  nozzle 

X  =  Y  -  (1 .0714  *  (Y  **  2)) 
1.5/(ratio  of  specific  heats)  =  1.0704 

-  compute  air  mass  flow  rate,  this  equation  is 
accurate  as  long  as  Y<0. 1 

MAIR  =  8.02*CD*PA*AREA*SQRT(X/(53.34*(459.7  +  TAIR))) 

AREA  =  area  of  air  inlet  nozzle 

TAIR  =  stagnation  temperature  of  air 

53.34  =  R  =  gas  constant 

10.  Air-to-Fuel  Ratio  (AF) ;  air  mass  flow  rate  divided  by 
fuel  mass  flow  rate 

AF  =  MAIR/MFUEL 
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11.   Scavenge  Ratio  (SCAV); 

SCAV  =  (MAIR  *  685  *  (TIN  +  459))/(PEXH  *  RPM) 
PEXH  =  exhaust  gas  pressure 
TIN  =  inlet  air  temperature 
685  =  (R/VD)  *  720 
VD  =  displacement  volume  =  53.0685  in 
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